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Abstract 

We explain the emergence and robustness of intense jets in highly turbulent plan- 
etary atmospheres, like on Jupiter, by a general approach of statistical mechanics of 
potential vorticity patches. The idea is that potential vorticity mixing leads to the 
formation of a steady organized coarse grained flow, corresponding to the statistical 
equilibrium state. Our starting point is the quasi-geostrophic 1-1/2 layer model, and 
we consider the relevant limit of a small Rossby radius of deformation. Then narrow 
jets are obtained, scaling like the Rossby radius of deformation. These jets can be 
either zonal, or closed into a ring bounding a vortex. Taking into account the effect of 
the beta effect and a sublayer deep shear flow, we predict an organization of the turbu- 
lent atmospheric layer into an oval-shaped vortex amidst a background shear. Such an 
isolated vortex is centered over an extremum of the equivalent topography (determined 
by the deep shear flow and beta-effect). This prediction is in agreement with analysis 
of wind data in major Jovian vortices (Great Red Spot and Oval BC). 

1 Introduction 

Atmospheric and oceanic flows are often organized into narrow jets. They can zonally flow 
around the planet like the jet streams in the earth stratosphere, or the eastward jet at 24 
in the northern hemisphere of Jupiter ( Maxworthy 1984). Jets can alternatively organize 
into rings, forming vortices, like the rings shed by the meandering of the Gulf-Stream in 
the western Atlantic Ocean. The flow field in Jupiter most famous feature, the Great 
Red Spot, is an oval-shaped jet, rotating in the anticyclonic direction and surrounding an 
interior area with a weak mean flow ( Dowling and Ingersoll 1989), see figure 1(a). Robust 
cyclonic vortices have a similar jet structure ( Hatzes et al 1981 ), see figure 1(b). 

Such jets and vortices are in a turbulent surrounding, and the persistence of their 
strength and concentration in the presence of eddy mixing is intriguing. The explanation 
proposed in this paper is based on a statistical mechanical approach: the narrow jet or 
vortex appears as the most probable state of the flow after a turbulent mixing of potential 
vorticity, taking into account constraints due to the quantities conserved by the dynamics, 
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especially energy. Such a statistical theory has been first proposed for the two-dimensional 
Euler equations by Kuz'min (1982), Robert (1990) , Robert and Sommeria (1991), Miller 
(1990). See Brandt et al (1999) for a recent review and discussion. This theory predicts an 
organization of two-dimensional turbulence into a steady flow (with fine scale, 'microscopic' 
vorticity fluctuations). Complete vorticity mixing is prevented by the conservation of the 
energy, which can be expressed as a constraint in the accessible vorticity fields. A similar, 
but quantitatively different, organization had been previously obtained with statistical me- 
chanics of singular point vortices with the mean field approximation, instead of continuous 
vorticity fields (Onsager 1949, Joyce and Montgomery 1973). 

Extension to the quasi-geostrophic (QG) model has been discussed by Sommeria et 
al (1991), Michel & Robert (1994), Kazantsev Sommeria and Verron (1998). This model 
describes a shallow water system with a weak vorticity in comparison with the planetary 
vorticity (small Rossby number), such that the pressure is in geostrophic balance, and the 
corresponding free surface deformation is supposed small in comparison with the layer thick- 
ness. For Jupiter the free surface would be rather at the bottom of the active atmospheric 
layer, floating on a denser fluid, as discussed by Dowling and Ingersoll (1989), see Dowling 
(1995) for a review. The gradient of planetary vorticity is accounted by a beta-effect. An 
additional beta-effect, depending on the latitude coordinate y, is introduced to represent 
the influence on the active atmospheric layer of a steady zonal flow in the deep interior, as 
discussed by Dowling and Ingersoll (1989). 

The free surface deformability, representing the strength of the density stratification, is 
controlled by the Rossby radius of deformation R*. The two-dimensional Euler equation is 
recovered in the limit of very strong stratification for which R* — > oo. We shall consider 
in this paper the opposite limit of weak stratification for which R* is much smaller than 
the scale of the system L. This limit is appropriate for large scale oceanic currents, as the 
radius of deformation is typically 10-100 km. For Jupiter, R* is estimated to be in the 
range 500-2500 km, while the Great Red Spot extends over 20,000 km in longitude, and 
10,000 km in latitude, so the limit R*/L — > seems relevant. We show that in this limit 
the statistical equilibrium is made of quiescent zones with well mixed uniform potential 
vorticity, bounded by jets with thickness of order R* . This provides therefore a general 
justification of jet persistence. Some of the ideas used have been already sketched in 
Sommeria et al (1991), but we here provide a systematic derivation and thorough analysis. 

In principle, the Quasi Geostrophic approximation breaks down for scales much larger 
than the radius of deformation, so that the limit R* /L — > seems inconsistent with the QG 
approximation. However the relevant scale is the jet width, which remains of order R*, so 
that the Quasi Geostrophic approximation remains valid in this limit. This point has been 
discussed by Marcus (1993) for the Great Red Spot, which he supposes to be a uniform 
Potential Vorticity ( PV ) spot surrounded by a uniform Potential Vorticity background ( 
we here justify this structure as the result of Potential Vorticity mixing with constraints on 
the conserved quantities ). Analyzing wind data in the Great Red Spot, Dowling & Ingersoll 
(1989) concluded that the QG approximation is good up within typically 30% error, which is 
reasonable to a first approximation. Statistical mechanics of the more general shallow water 
system (to be published), predicts a similar jet structure. The present Quasi Geostrophic 
results therefore provide a good description as a first approximation. 

We first consider the case without beta-effect in section 2. We furthermore assume 
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periodic boundary conditions (along both coordinates) in this section to avoid consideration 
of boundary effects. Starting from some initial condition with patches of uniform PV, we 
find that these patches mix with uniform density (probability) in two sub-domains, with 
strong density gradient at the interface, corresponding to a free jet. The coexistence of the 
two sub-domains can be interpreted as an equilibrium between two thermodynamic phases. 
We find that the interface has a free energy per unit of length, and its minimization leads to 
a minimum length at equilibrium. This results in a constant radius of curvature, in analogy 
with surface tension effects in thermodynamics, leading to spherical bubbles or droplets. 
The range of the vortex interaction is of the order R* , therefore becoming very small in 
the limit of small radius of deformation, so the statistical equilibrium indeed behaves like 
in usual thermodynamics with short range molecular interactions. This contrasts with 
the case of Euler equation, with long range vortex interactions, analogous to gravitational 
effects (Chavanis Sommeria and Robert 1996, Chavanis 1998). 

Figure |?] summarizes the calculated equilibrium states, depending on the total energy 
and a parameter B representing an asymmetry between the initial Potential Vorticity patch 
areas, before the mixing process. We obtain straight jets for a weak asymmetry and circular 
jets for higher asymmetry. Such circular a jet reduces to an axisymmetric vortex, with 
radius of order R* , in the limit of low energy. 

We discuss the influence of the beta-effect or the deep zonal flow in section 3. The 
channel geometry, representing a zonal band periodic in the longitude x is appropriate for 
that study. With the usual beta-effect (3y, linear in the transverse coordinate y, statistical 
equilibrium is, depending on the initial parameters, a zonal flow, or a meandering eastward 
jet, or a uniform velocity v m = R 2 j3 whose induced free surface slope cancels the beta-effect 
( uniformization of Potential Vorticity ) on which circular vortices can coexist. 

For more general beta-effects, due to the deep zonal flow, we find that the jet curvature 
depends on latitude y. In particular a quadratic beta effect ay 2 leads to oval-shape jets, 
similar to the Great Red Spot. Using the determination of the sublayer flow from Voyager 
data by Dowling and Ingersoll (1989), we show in section 4, that such a quadratic effective 
beta-effect is indeed a realistic model for Jupiter atmosphere in the latitude range of the 
Great Red Spot and the White Ovals, the other major coherent vortices on Jupiter. Using 
these data on beta-effect, as well as the shear in the zonal flow at the latitude of the Great 
Red Spot, the jet width and its maximum velocity, we deduce all the parameters of our 
model. 

2 The case with periodic boundary conditions 

2.1 The dynamical system 

We start from the barotropic Quasi Geostrophic (QG) equation : 



Oq 

di 



+ v • Vq = 



(1) 
(2) 
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v = — z A Vip 



(3) 



where q is the potential vorticity ( PV ), advected by the non-divergent velocity v, i/j is 
the stream function, R is the internal Rossby deformation radius between the layer of fluid 
under consideration and a deep layer unaffected by the dynamics, x and y are respectively 
the zonal and meridional coordinates ( x is directed eastward and y poleward ). The 
term h(y) represents the combined effect of the planetary vorticity gradient and of a given 
stationary zonal flow in the deep layer, with stream function ipd(y) : h(y) = —(3y + ipd/R 2 . 
This deep flow induces a constant deformation of the free surface, acting like a topography 
on the active layer. We shall therefore call h(y) the 'topography', and study its influence 
in section 3. Let us assume h(y) = in this section. We define the QG equations dlL@) 
in the non-dimensional square D = [— |, ^] 2 . R is then the ratio of the internal Rossby 
deformation radius R* to the physical scale of the domain L. 

Let (/) = f D fd 2 r be the average of / on D for any function /. Physically, as the 
stream function i/j is related to the geostrophic pressure, (if)) is proportional to the mean 
height at the interface between the fluid layer and the bottom layer, and due to the mass 
conservation it must be constant (Pedlosky 1987). We make the choice 

m = o (4) 

without loss of generality. 

The total circulation is (q) = (—Aijj + i^/R 2 ) = (tp/R 2 ) due to the periodic boundary 
conditions. Therefore 



(q) = (5) 

We note that the Dirichlet problem (|2|) on D with periodic boundary conditions has a 
unique solution ip for a given PV field. 

Due to the periodic conditions for tp, the linear momentum is also equal to 0, 



The energy 



E = - 
2 



(v) = (6) 



\ j D q^d 2 v = \j D [m) 2 + f [2 ]d 2 r (7) 



is conserved ( we note that the first term in the right hand side of (0) is the kinetic energy 
whereas the second one is the gravitational available potential energy ). 
The integrals 

Cf(q) = f f(q)d 2 v (8) 

for any continuous function / are also conserved, in particular the different moments of the 
PV. In the case of an initial condition made of a finite number of PV levels, the areas initially 
occupied by each of these levels is conserved, and this is equivalent to the conservation of 
all the constants of motion M) . 
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2.2 The statistical mechanics on a two PV levels configuration. 



2.2.1 The macroscopic description. 

The QG equations (|l|) (|2|) are known to develop very complex vorticity filaments. Because 
of the rapidly increasing amount of information it would require, as time goes on, a de- 
terministic description of the flow for long time is both unrealistic and meaningless. The 
statistical theory adopts a probabilistic description for the vorticity field. The statistical 
equilibrium depends on the energy and of the global probability distribution of PV levels. 
Various previous studies (Sommeria & al 1991), (Kazantsev & al 1998) indicate that a 
model with only two PV levels provides a good approximation in many cases. The determi- 
nation of the statistical equilibrium is then simplified as it depends only on the energy, on 
the two PV levels, denoted q = a\ and q = a_i and on their respective areas A and (1 — A) 
in D. The number of free parameters can be further reduced by appropriate scaling. Indeed 
a change in the time unit permits to define the PV levels up to a multiplicative constant, 
and we choose for the sake of simplicity : 

2^-1 (9) 
and define the non-dimensional parameter B as : 

(io) 

The condition (||) of zero mean PV imposes that a\A + a^\(l — A) = 0. This means that a\ 
and a_i must be of opposite sign and, using (|) and flu]), A = (1 - B)/2. The distribution 
of PV levels is therefore fully characterized by the single asymmetry parameter B, which 
takes values between -1 and +1. The symmetric case of two PV patches with equal area 
A = 1/2 corresponds to B = 0, while the case of a patch with small area (but high PV, 
such that ( q ) = 0) corresponds to B — ► 1. Note that we can restrict the discussion to 
B > 1 as the QG system is symmetric by a change of sign of the PV. 

The two PV levels mix due to turbulent effects, and the resulting state is locally de- 
scribed by the local probability (local area proportion) p(r) to find the first level at the 
location r. The probability to find the complementary PV level a_i is 1— p, and the locally 
averaged PV at each point is then 

q(r) = a lP (r) + o_i(l - p(r)) =2^-^+5 (11) 

where the second relation is obtained by using d) and (0). 

As we consider the evolution of two PV patches, the conservation of all invariants (^) 
is equivalent to the conservation of the area A of the patch with PV value a\ ( the area of 
the other PV level a_i being 1 — A). The integral of p over the domain must be therefore 
equal to the initial area A ( the patch with PV level a\ is mixed but globally conserved ), 

A = ±=JL = [ p(r)d 2 r (12) 
* Jd 
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As the effect of local PV fluctuations is filtered out by integration, the stream function 
and the velocity field are fully determined by the locally averaged PV q as the solution of 

ib 

q = —Aip + — r ; ip periodic (13) 
R z 

and v = — z A V?/> 
Therefore the energy is also expressed in terms of the field q : 



2 



D 



2 



d 2 r = \ I q^d 2 r (14) 



2 



D 



Here the energy of the 'microscopic' PV fluctuations has been neglected (replacing q by 
q), as justified in the case of Euler equation by Robert and Sommeria (1991). Indeed, 
considering a 'cutoff' for the microscopic fluctuations much smaller than R, the small scale 
dynamics coincides with the Euler case. 

The central result of the statistical mechanics of the QG equations (Qj2|) is that, under 
an ergodic hypothesis, we expect the long time dynamics to converge towards the Gibbs 
states defined by maximizing the mixing entropy 

S = - / [ p(r) lnp(r) + (1 - p(r)) ln(l - p(r)) }d 2 r (15) 
Jd 

under the constraints of the global PV distribution (12) and energy (|14|). It can be shown 
that the microscopic states satisfying the constraints given by the conservation laws are 
overwhelmly concentrated near the Gibbs state, which is therefore likely to be reached 
after a complex flow evolution. A good justification of this statement is obtained by the 
construction of converging sequences of approximations of the QG equation (|l|j2|), in finite 
dimensional vector spaces, for which a Liouville theorem holds. This is a straightforward 
translation of the work of Robert (1999) for 2D Euler equations. The sequence of such 
Liouville measures has then the desired concentration properties as (0j2|) enters in the 
context considered in Michel & Robert (1994). 

2.2.2 The Gibbs states 

Following Robert & Sommeria (1991), we seek maxima of the entropy ( |i~5|) under the 



constraints (12) and Ql4|) . To account for these constraints, we introduce two corresponding 
Lagrange multipliers, which we denote 2a and —C/R 2 for convenience in future calculations. 
Then the first variation of the functionals satisfies : 

SS - 2aSA + %r5E = 
R z 

for all variations 5p of the probability field p. After straightforward differentiation we 
obtain: 

5S = - [ [ lnp - ln(l - p) ]5pd 2 r , 5 A = [ 5pd 2 v 

JD JD 

6E= [ i)5qd 2 r = / 2ip5pd 2 r (16) 

JD JD 
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where the expression of 5E has been obtained by integrating by part and expressing q by 
(JTT|) . Then we can write the first variation under the form J D [ — lnp + ln(l — p) — 2a + 
2Ci/)/R 2 ] dp d 2 r which must vanish for any small variation Sp. This implies that the 
integrand must vanish, and yields the equation for the optimum state: 

1 - tanhfo - S£) 
P= \ (17) 

and using ( pT| ) and ([T3|), the partial differential equation 

q = -Aip + ^ = B - tanh (a - °*\ (18) 

determining the Gibbs states (statistical equilibrium). Prom now on we forget the q over- 
line for the locally averaged PV and refer to it as the PV. 

Therefore, we have shown that for any solution of the variational problem, two constants 
a and C exist such that ip satisfies ( JT8|) . Conversely it can be proved that for any such 
two constants, a solution to equation (jig), in general not unique, always exists. Then p 
associated with one of these solutions by (|l7]) is a critical point of the 'free energy' — S(p) + 
2aA(p) — -UiE{p) (i.e. its first variation vanishes). Then the Lagrange multipliers are not 



given but have to be calculated by prescribing the constraints ( |12; ) and (14) corresponding 
to the two parameters B and E respectively, given by the initial condition. Furthermore, 
among the states of given energy E and asymmetry parameter B, we shall select the actual 
maxima. 

Finally, let us find a lower bound for the parameter C of the Gibbs states with non-zero 

energy (i.e. if) is not constant on D). Multiplying (Il8|) by —Aip, integrating by part and 

(h. 
K 2 



defining f(Ctp) = B — tanh(o — 7^), we obtain : 



C 



J ((Aip) 2 + ^(VVQ 2 ) cPr 



f D -f'(CiP)(ViP) 2 d 2 r 
From which, using < —f'(Cip) < it follows that when ip is not constant : 

C > 1 (19) 



2.3 The limit of small Rossby deformation radius 

As suggested by oceanographic or Jovian parameters, we seek solutions for the Gibbs states 
equation in the limit of a small ratio between the Rossby deformation radius and the length 
scale of the domain : R << 1 with our non-dimensional coordinates Q 

1 Modica (1987) considered the minimization of the functional E e (u) = J n [e (Vu) 2 + Wo(«)]dx with the 
constraint J n «(x)dx = m in the limit e — > + where Wo is a real function with two relative minima. He 
proved, in a mathematical framework, working with bounded variations functions, that if (« £ ) are solutions 
of this variational problem, for any subsequence of (u e ) converging in _L 1 (fi) as e — > 0, this subsequence 
converge to a function uo which takes only the values where Wo reaches its minima ; with the interface 
between the corresponding subdomains having minimal area ( See Modica (1987) for a precise statement ). 

We note that the Euler equation of this variational problem may be the same as the Gibbs States equation 
( H ) for a convenient choice of Wo- However as the variational problem itself is different this beautiful result 
cannot be used in our context. 
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2.3.1 The uniform subdomains 



Then we expect that the Laplacian in the Gibbs states equation (|T^) can be neglected with 
respect to ip/R 2 , except possibly in transition regions of small area. This transforms jl8| ) 
into the algebraic equation : 



Q = ^2=B - tanh (a - ^ (20) 

Depending on the parameters, this equation has either one, two or three solutions, 
denoted f/'-ij^'o an d ipi in increasing order (see figure || ). The case with a single solution 
would correspond to a uniform ip, which should be equal to due to the condition (ip) = 0. 
This is only possible for E = 0. Otherwise, we have therefore two or three solutions, with 
different solutions occurring in subdomains. This condition of multiple solutions requires 
that the maximum slope for the right hand side of ( p0| ) must be greater than 1/R 2 ; this is 
always realized due to the inequality (|T^). Furthermore a must be in an interval centered 
in CB ( a = CB in the symmetric case of figure ^ ). 

At the interface between two constant stream function subdomains, a strong gradient of 
ip necessarily occurs, corresponding to a jet along the interface. These jets give first order 
contributions to the entropy and energy, but let us first describe the zero order problem. 
Suppose that ip takes the value ip\ ( resp ip-\ ) in subdomains of total area A\ ( resp 
A-i ). The reason why we do not select the value ipo will soon become clear. Using 
(|Tl| ) we conclude that the probability p takes two constants values p±\ in their respective 
subdomains. The two areas A±i ( measured from the middle of the jet ) are complementary 
such that : 



Ai + A_ x = l (21) 
Furthermore the constraint (|j) of zero domain average for ip implies at zero order, 

V^+V-iA.! =0, (22) 
or equivalently, using q±i = i{j±i/R 2 , (|Tl|) and (pi]) : 

2A l L 2(1 - A x ) (p-x ~l)=-B (23) 



This can be obtained as well from the constraint on the (microscopic) PV patch area ([12]). 
The energy inside the subdomains reduces to the potential term ip 2 /2R 2 , since velocity 
vanishes. This area energy Ea can be computed in terms of p±\ using q±\ = ip±i/R 2 and 

Ea - ^ lAl + 2 £ lA - 1 = AMpi) + (1 - AMP-i) 

( ( 1\ 2 / 1\ B 2 \ (24) 

with e (p)=R 2 [2ip--\ +2B[p--\+—\ 
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There is also an energy in the jet at the interface of subdomains, but it is small with respect 
to E A . Indeed the velocity in the jet, of width R, is of order (ifti — ip_\)/R ~ ip/R, and the 
corresponding integrated kinetic energy is of order tp 2 /R. This is small in comparison with 
the area energy Ea (mostly potential) of order ip 2 /R 2 . A precise calculation will confirm 
this estimate in next sub-section. 

We need to determine three unknown, the area Ai and the probabilities p±i of the PV 
level a\ in each subdomain, while the constraints (^) and (|24] ) provide two relations. An 
additional relation will be given by entropy maximization. As we neglect the jet area, the 
entropy reduces at order zero to the area entropy : 

S A = Axsipt) + (1 - Ai)s(p_i) 
with s{p) = —plogp — (1 — p) log(l — p) 



Thus the zero order problem corresponds to the maximization of the area entropy (25) 
with respect to the 3 parameters p±\ and A\, under the 2 constraints ( p3[ ) and (p4j). A 
necessary condition for a solution of this variational problem is the existence of two Lagrange 
parameters ao and Co, associated respectively with the circulation constraint (^) and with 
the energy constraint d2J), such that the first variations of the total free energy, 

F A = -S A -^E A + a ^-, (26) 



vanish. Let us calculate F A using ( |23| ) and (f24|): 
F A = A 1 f(p 1 ) + (l-A 1 )f(p- 1 ), 
with /(p) = -s(p) - 2C (p-^j " 2 (C B - a ) (p - ^ - ^ - a Q B. 2? 

The vanishing of the variations with respect to p\ and p_i gives that f(p±i) are local 
minima of the free energy f(p). It is easily proven that the function / has two local minima 
and one local maximum ( for Co > 1 and (CqB — ao) small enough ) ( see figure ||). The local 
maximum is achieved for po corresponding to the value i/jq . It is the reason why it has not 
been taken into account in this analysis. In addition, the vanishing of the first variations 
with respect to the area A% imposes the free energies f(p±i) in the two subdomains to 
be equal. This is like the condition of thermodynamic equilibrium for a chemical species 
shared by two coexisting phases. 



In the expression (27) of f(p), the entropy term s(p) is symmetric with respect to 
p = |, as well as the quadratic term. Therefore if the linear term in (p — |) vanishes the 
two maxima are equal, with p±i symmetric with respect to \. The addition of a linear 
term obviously breaks this condition of two equal maxima, so the coefficient of the linear 
term must vanish, thus : 

a = C B. (28) 
Since p±i are symmetric with respect to ^, we introduce the parameter u by : 

p± 1 = i(l±u). (29) 
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Using (|TT|) , (|23|) we deduce: 



V>±i = R 2 (B±u) 



(30) 



From ( p2[ ) we state that the two constant stream function ( p0| ) have to be of opposite sign, 
so that u > |£?|. Introducing (H) in the circulation constraint (p3|), and using (|2l|) , we get 

A ±1 = 5 (l =F f ) ■ (3D 
Using these results, the energy ( |24| ) becomes 

E^E A = ^-(u 2 -B 2 ) (32) 

This relates the parameter u to the given energy E and asymmetry parameter i?. Finally 
the condition that f(p±i) are maxima of / leads to : 

u = tanh(Cou), (33) 

which determines the 'temperature' parameter Co, as represented in figure |j. Therefore all 
the quantities are determined from the asymmetry parameter B and from the parameter 
u, related to the energy by (|32|) . 



In the limit of low energy, u — > \B\, when for instance B > 0, then A\ goes to zero, 
so that ij)—x tends to occupy the whole domain. This state is the most mixed one com- 
patible with the constraint of a given value of B (or equivalently a given initial patch area 
A = (1 — B)/2). In the opposite limit u — * 1, we see from ( |30"| ) that in the two subdomains 
q = ip/R 2 tends to the two initial PV levels a\ = 1 + B and a_i = —1 + B. Thus, this state 
is an unmixed state. It achieves the maximum possible energy E = \{l-B 2 ) under the 
constraint of a given patch area. We conclude that the parameter it, or the related 'tem- 
perature' Co, linked with the difference between the energy and the maximum accessible 
energy for the two given initial PV levels, characterizes the mixing of these two PV levels. 
We shall call u the segregation parameter, as it quantifies the segregation of the PV level 
a± ( or its complementary a_i ) between the two phases. 



Let us now study the interface between the subdomains. 



2.3.2 Interior Jets 

At the interface between two constant stream function subdomains, a strong gradient of ifi 
necessarily occurs, corresponding to a jet along the interface. To study these jets, we come 
back to the Gibbs state equation fll8|). We expect the Lagrange parameters a and C to be 
close to the zero order parameters «o and Co, computed in the previous sub-section, so we 
use a = ao and C = Co to calculate the jet structure. In such a jet, we cannot neglect the 
Laplacian term in (18), but a boundary layer type approximation can be used: we neglect 
the tangential derivative with respect to the derivative along the coordinate normal to the 
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interface, £. Accordingly, we neglect the inverse of the curvature radius of the jet with 
respect to 1/R. 

Thus, from the Gibbs states equation ([i8[), using (|28|), we deduce the jet equation : 



** + *-B-t**(c(B*]] (34) 



d( 2 R 2 V V R 2 

As the stream function depends only on the normal coordinate the velocity is tangent to 
the interface, forming a jet with a typical width scaling like R. We thus make the change 
of variables defined by : 

r=^i <P^-B + ^, (35) 
leading to the rescaled jet equation: 

^| = -tanh(C <j>) + <i> (36) 

The jet equation ([36]) is similar to a one dimensional equation of motion for a particle 
( with the position depending on a time r) under the action of a force —dU/dcj) deriving 
from the potential, 

= ln(cosh(C ^) _ ^ 
Co 2 



1 "- 2 



represented in figure |2|(b) . In its trajectory the particle energy is conserved : 

+ {/(</>) = Cst (38) 



2 \dr 

Let <pi = ifii/ R 2 — B, i = — 1, 0, 1, corresponding to the solutions tpi of the algebraic equation 



(|2C)|). From (|3C|), we have (j)±\ = ±u. Note that the stationary limit of (36), which must 
be reached for lim^-too, yields again (|3lf). Moreover, the particle energy conservation (35) 
imposes the integrability condition, 

U(^) = Ufa), (39) 
which is indeed satisfied due to the symmetry of the potential U. We note that the Lagrange 



parameter determination (28) and the symmetry of the probabilities (^) with respect to 
7j could have been deduced from this integrability condition (|39| ) instead of minimizing 



the free energy ( p7| ) (we shall proceed in this way in section 3 to take into account the 
beta-effect). 

The jet equation ( plf ) has been numerically integrated. Figure || shows a typical stream 
function and velocity profile in the jets. Figure ||(a) shows how the jet width depends on 
the segregation parameter u. We note that the width of the jet is an increasing function 
of the mixing and therefore a decreasing function of the energy. Figure |^(b) shows the 
dependence in u of the total non-dimensional energy e(u) = \ (dcp/dr) dr and of the 
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maximum non dimensional jet velocity (dcfr/dr) 



As the jet structure (36) does not depend on the coordinate tangent to the jet, we can 
define the jet entropy ( respectively energy, free energy ) per unit length S j e t ( respectively 
Ejet, Fjet )• Multiplied by the jet length, these quantities are the first order corrections 
to the entropy ( respectively energy, free energy ). Using the change of variables (|35|), we 
calculate the jet entropy per unit length : 



/+oo 
[s{p{t)) - s{p±i)]dT 
-oo 



where s is defined in (p5|), and p±\ are defined in (|29|), Using the probability equation (|T 
and (^) we obtain : 

/+oo 
[*(</>)) - ~s(4>±i))]dT (40) 
-oo 

involving the function s((p) = ln(cosh(Co0)) — Co</>tanh(Co</>). Similarly we straightfor- 
wardly calculate the potential and kinetic energy per unit length for the jet : 

" T L ( * - 4)dT E} " = T L (5F J rfT 

We use the integral p^ ) to calculate d4>/dr and conclude : 

/+oo 
[sr(0)-S(0±i)]<fc-> (41) 
-oo 

with g((/>) = — in(co^(Cb^) _|_ -Q ue |- S y mme t r y Q f ^jjg jets, the jets provide no 
perturbation to the zero-order circulation, so there is no circulation term in the jet free 
energy expression : F Jet = -Sj et - C /R 2 Ej et . Then 

/+oo 
H<t>) - K<t>±i)]dT (42) 
-oo 

with h(4>) = -<p{4> - tanh(Co^)). 

Let us study the sign of Fj et . As <pi verifies 4>\ = tanh(Co</>i) and as </>(r) is an increasing 
function of r with lim T ^ +oc 4>{t) = 4>\ we conclude that h{4>) — h{4>i) > for any r > 0. 
Thus Fj e t > 0. Using the analogy with usual thermodynamics, the 'surface tension' is 
positive. This favors large 'bubbles' which minimize the interfacial length and therefore 
the corresponding free energy (p7|). Our initial hypothesis of well separated domains with 
uniform ifi is thus supported, as discussed more precisely in next subsection. 

2.3.3 Selection of the sub-domain shape 

The above analysis has permitted us to determine the areas of sub domains on which the 
stream function tp takes the constant values ip±i, but the subdomains shape is still to be 
selected. There is an analogy with two phases coexisting in thermodynamic equilibrium, for 
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instance a gas bubble in a liquid medium, for which a classical thermodynamic argument 
explains the spherical shape of the bubble by minimizing its free energy, proportional to 
the bubble area. Our system is isolated rather than in a thermal bath, but the jet energy 
is small (of order R) with respect to the total energy. Therefore the subdomain interior 
behaves like a thermal bath with respect to the jet, so the usual argument on free energy 
minimization applies. We shall now show this more precisely by directly maximizing the 
total entropy with constraints, taking into account the jet contribution. 

The jet with length L has an entropy Sj e tL and energy Ej et L. Since the total energy 
E = E A (C) + LEj e t is given, the jet has also an indirect influence in the area energy 
E A . This small energy change 5E A results in a corresponding change in the area entropy 
5S A = -{C /R 2 )5E A , from the condition (gjgg ) of zero first variations. Note that there is 



no area change 5 A since the jet is symmetric and has therefore no influence in the condition 
( |l2| ) of a given integral of p (the difference in p with the case of two uniform patches with 
boundary at the jet center has zero integral). Therefore, adding the direct and indirect 
contribution of the jet entropy leads to the total entropy 

S = S A (C ) + (j^Ejet + S Je ^j L = S A (C ) - F Jet L (43) 

where <SU(Co) is the zero order, area entropy, obtained in the limit of vanishing jet width ^] 



We deduce from (42) that the maximization of the entropy is achieved by minimizing 
the total free energy Fj e tL, which we have proved to be positive at the end of previous 
sub-section. Thus we conclude that the maximum entropy state minimizes the jet length, 
with a given area of the subdomains (|3l|). The subdomains shape will therefore be a circle 
or a stripe. More precisely if A\ < 1/tt the jet forms a circle enclosing the positive constant 
stream function domain ( the jet bounds a cyclonic vortex ), if 1/tt < A\ < 1 — 1/tt two 
straight lines jets form a stripe and if A% > 1 — 1/tt the jet form a circle enclosing the 
negative constant stream function domain ( the jet bounds anti-cyclonic vortex ). 

The different types of solutions can be summarized in a (E,B) diagram : figure [7[ The 
outer parabola is the maximum energy achievable for a fixed B : E = R 2 (l - B 2 )/2. The 
frontier lines between the straight jets and the circular jets corresponds to A\ = 1/tt or 
A_i = 1/tt. It has been calculated using (^TJ) and @) : E = R 2 B 2 (2tt - 2)/(vr - 2) 2 . Note 
that the maximum accessible energy is in R 2 , but it has been scaled by the normalization 
condition (]9|) on PV levels, so the real energy is not bounded. 

All this analysis assumes that the vortex size is much larger than the jet width A, given 



by figure y. In other words, the area Ai or A-\ (31) must be larger than (2A) . This is 
satisfied on the right side of the dashed line represented in figure 0. The dashed line itself 
corresponds to the equality, and the condition of large vortex is clearly not satisfied on its 
left side, for low energy. The position of the dashed line depends on the numerical value of 
R (it has been here numerically computed for R = 0.03), and it gets closer to the origin as 
R —* 0. We shall now determine the statistical equilibrium in this case of low energy. 

2 This reasoning to obtain the first order entropy can be precised by evaluating explicitly the first order 
modification of the Lagrange parameter C ( let say Ci = C — Co ) due to the jet energy, and the first order 
modification of the Lagrange parameter a ( let say a 1 = a — Qo ) due to the jet curvature and computing 
the first order entropy from its definition (Ifq). We have calculated the first order entropy in this way to 
actually obtain mM. 
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2.3.4 Axisymmetric vortices 



We have noted in subsection ( 2.3.1| ) that, when for instance B > 0, in the limit of small 



energy (E — > or equivalently u — ► \B\, for fixed B and R), the area A_i occupied by ip-i 
tends to 1, the whole domain. Therefore, in this limit, the complementary area A\ tends to 
and the vortex becomes smaller than the deformation radius, so we can no more neglect 
the curvature radius of the jet. 

In this limit u —* \B\, as the vortex has a small area with respect to the total domain, 
it is not affected by the boundary conditions, so it can be supposed axisymmetric. From 



the general Gibbs states equation (|18|), we deduce the axisymmetric vortex equation : 

-de~cdt + & = B " tanl T" «v m 

As R will be a typical scale of the vortex, we make the change of variable, 

C = Rr ; <j> = ~ + (45) 
leading to the rescaled axisymmetric vortex equation : 

-f*-±^ + * + £=S + tanh(C<fr)) (46) 
dr A r ar L 

From now on, we shall consider the case B > (the case B < is just the symmetric case 
of a negative vortex). 

For this equation to describe a localized vortex, we impose lim r -, too (f){r) = <j>—\ = 
—a/C + ip-i/R 2 , where ip-± is the positive solution of the algebraic equation Since 
nearly the whole fluid domain is covered by the asymptotic stream function "0_x outside 
the vortex, the condition of zero total circulation (q) = (ip)/R 2 = imposes that V-i — 
(it is of order R), so that 0_i = —a/C, and the algebraic equation ( |20|) then leads to: 

a = argtanh(£>) (47) 

We can thus eliminate a in (|4"6|), leading to an equation depending on two parameters, 
B and C, 

^ = -^-tanh ( ^) + 0-i? + ^^ 
dr £ r ar L 

^-{r = 0) = and lim <p(r) = — & 

dr r^oo C 

where the regularity condition at r = has been included. Let us consider, as in section 



( 2.3.2 ), the analogy of equation (fig) with a one particle motion with 'position' cj) and 'time' 



r. 



The last four terms on the right-hand side of ([46|) can be written as the derivative 
-dU/d(j) of the potential, 

m = ln(cosh(^)) _ £ / _ argtanh^ \ ^ (4g) 
O Z \ O / 
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(represented in figure ||), while the first term can be interpreted as a friction effect. Indeed, 
an integration of (|48|) leads to: 



-oo n / j ± \ 2 



- I7(0(r = 0)) = - Jf ^(j£J dr<0 ( 5 °) 

Thus, in figure ||(a), the hatched area on the right side must be greater than the one on 
the left (since (U((f>i) - f7(0_i)) > U((f>(r = 0)-U {<f>-\) > 0). It is clear from figure | that 
this is possible only if </> < and a/C < B, or, using (p7|), C > a/B = argtanhi?/S. 
The value C = a/B corresponds to the integrability condition ( |39|) when the effect of jet 
curvature is neglected. This effect is now taken into account by the departure of C from 
this value, which we shall denote AC = C — argtanhi?/£>. Then AC > and we expect to 
recover the results of section 2.3.2 in the limit AC —* 0. Moreover, we must reach a uniform 
stream function at large distance, solution of the algebraic equation (|20|), so it must have 
three solutions. We see in figure || that the corresponding AC must not exceed a maximal 
value, denoted AC max . 

We can prove that for any B > and argtanh£?/£> < C < argtanh .B/.B + AC max , 
equation (E3) has a unique solution. Such solutions have been numerically obtained for 
B = 0.75 and < AC < AC max . Corresponding stream function profiles are shown in 
figure ||. 

As AC is decreased from AC max to zero, two stages can be seen in figure [| First the 
maximum value for the stream function is increased while the mean width of the vortex 
remains of the order of R. In a second stage, when AC goes to zero, as we are closer to 



the integrability condition for big vortices (3£), (j) remains longer in the vicinity of <pi so 



the vortex size increases. Note that the energy monotically increases as AC is decreased, 
first by an increase in the vortex maximum stream function and then by an increase in size. 
Finally the case of a jet with negligible curvature studied in subsection 2.3.2 is reached 
when AC -> 0. 

In conclusion, we have shown that in the limit of small energy, with fixed B and R, the 
Gibbs states are approximated by axisymmetric vortices, whose radial structure depends on 
the parameter AC, which monotonically decreases from AC max to as energy is increased. 



2.3.5 The linear approximation for the Gibbs states 

The previous discussion of axisymmetric vortices was concerned with the limit of small 
energy with fixed B and (small) fixed R. We consider now the limit of small E and B, 
i.e. the neighborhood of the origin in the phase diagram of figure [?]. Then from fl32]), 
u — ► \B\ << 1. Figure || shows that for \u\ « 1, the jet width diverges and therefore the 
jet tends to develop on the scale of the whole domain, so the approximation of a localized 
jet, or an isolated axisymmetric vortex, falls down. 

In this limit of small E and B, we can however linearize the Gibbs states equation fll8|) , 
following the work of Chavanis and Sommeria (1996) for Euler equation. After linearization, 
solutions are expressed in terms of the eigenmodes of the Laplacian, and only the first 
eigenmodes can be entropy maxima. 

These results are unchanged by the linear deformation term ip / R 2 , so the work of Cha- 
vanis and Sommeria (1996) directly applies here. With the periodic boundary conditions, 
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the first eigenmode of the Laplacian, a sine function of one of the coordinates, for instance 
y, is thus selected. This corresponds to the low energy limit of the two jet configuration 
shown in figure [?]. The next eigenmode, in sin(-7rx) sin(Try), has the topology of the vortex 
states. A competition between these two modes is expected in the neighborhood of the 
origin for small E and B. Note that the range of validity of the linear approximation is 
limited to a smaller range of parameters than in the Euler case, and this range of validity 
gets smaller and smaller as R — > 0. The dominant solution with uniform subdomains and 
interfacial jets relies by contrast on the tanh like relation between PV and stream function, 
and it is genuinely non-linear. 

3 The channel case 

We now consider the channel geometry, which represents a zonal band around a given 
latitude. It is then natural to introduce a beta-effect, or a mean sublayer zonal flow ( 
topography ), with the term h(y) in (|l|). We shall study the two cases of a linear h(y) 
( beta effect and/or uniform velocity for the sublayer flow ) or a quadratic h(y). We 
shall follow the presentation for the periodic boundary conditions, stressing only the new 
features. 

3.1 The dynamical system 

Let us consider the barotropic QG equations ([l], p] and ^) in a channel D = \—\, \] 2 with 
the velocity v tangent to the boundary for y = ±2 and 1-periodicity in the zonal direction. 
Thus we choose for the boundary conditions a constant ip, denoted tpb, the same on the two 
boundaries y = We note that, due to these conditions, the physical momentum (J|) is 
equal to zero. It is always possible to satisfy this condition by a change of reference frame 
with a zonal velocity V such that it moves with the center of mass of the fluid layer, and a 
corresponding change of the deep flow, resulting in an additional beta effect h — ► h + 

As in section (^), we need to specify the gage constant in the stream function ip, and 
we generalize the integral condition (Q) as, 

j^r - <%)> = 0. (51) 

The total mass (ip) is then constant in time (but not the boundary value ipb in general). 
With these conditions, the Dirichlet problem (||) has a unique solution if) for a given PV 
field q. We note that the scale unit is chosen such that the area of D is equal to 1. 

The integral of any functions of the potential vorticity (||) is still conserved. Let in 
particular T be the global PV, or circulation : 

T=(q)= [ -Aij) d 2 r = [ v.cfl (52) 
J D J 3D 

By contrast with the doubly periodic boundary conditions, the circulation T is not neces- 
sarily equal to zero. The expression of the energy in terms of the PV (see equation ([?])) is 
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therefore modified (due to the boundary term in the integration by parts) : 

d 2 r = ^J^ q + h(y))^d 2 r- l -T^ b (53) 

Due to the invariance under zonal translation of the system, another conserved quantity 
exists : 

M = yqd 2 r (54) 
Jd 

This constant moment fixes the 'center of mass' latitude for the PV field. 
3.2 General form of the Gibbs states 

et us consider the statistical mechanics on a two PV level configuration : the initial states 
is made of patches with two levels of potential vorticity, q = a\ and q = a_i, occupying 
respectively the areas A and (1 — A) in D. We keep the normalization (|9|) and definition 
(|i~0| ) for B. Now, since the circulation V is non-zero, the area A is related to B by A = 
(1 — B)/2 + r/2. The boundary term in the expression of the energy ( |53|) leads to an 
obvious change in the energy variation (|l^). Let 7 be the Lagrange multiplier associated 
with the conservation (^4|) of the momentum M. 

Adapting the periodic case computations, we then calculate the probability equation 
and the Gibbs state equation : 

1 - tanh ( a' - ^ + jy) 
P = L (55) 



E = \ 
2 .id 



q = -Atp + Ky) = B — tanh (a' - ^ + yy) (56) 

with a' = a + Cipb/ R 2 ■ These results generalize ( |l7l) and (|l8l) of section 2. 

In the case of a Gibbs state depending on x, the Lagrange parameter 7 is related to a 
zonal propagation of the equilibrium structure. The statistical theory only predicts a set 
of equilibria shifted in x, but introducing the result back in the dynamical equation yields 
the propagation. Indeed the Gibbs state equation ([56]) is of the form q = f(ip,y), which 
can be inverted (as it is monotonous in ifj, see Robert and Sommeria (1991)) to yield: 

^ = 9{q)+R 2l -§ (57) 

where g is a function of the potential vorticity. From this relation we calculate the velocity 
using (||) : v = i? 2 ^x — g'(q)z A Vq. As PV is advected ( equation ([!]) ) we obtain : 

Thus the PV field is invariant in a frame propagating with the zonal speed V sr = R 2 ^. 
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3.3 The limit of small Rossby deformation radius. 



In this sub-section we propose to analyze the Gibbs state equation (|56|) in the limit of small 
deformation radius ( R « 1 ). The main difference with the periodic case resides in the 
latudinally depending topography h(y), resulting in two effects. Firstly the subdomains of 
uniform PV are no more strictly uniform, and contain a weak zonal flow. Secondly the 
jet curvature is no longer constant in general, but depends on the local topography. Using 
the same boundary layer approximation as in the periodic case, the Laplacian term in the 
Gibbs state equation ( |5"6| ) will be neglected, except possibly in an interior jet and in the 
vicinity of the boundaries y = ±5 ( boundary jets ). 

Outside such jets, ([56]) reduces to an algebraic equation, like previously : 



- Ky) = B - tanh (a' ~C^+ 7 y) • (59) 

This is like (p0[), replacing the constants B by B + h{y) and a by a' + ^y. The three 
solutions can be still visualized by figure ^, but the position of the straight line with respect 
to the tanh curve now depends on y, due to the terms h{y) and 7?/. We assume that this 
dependence is linear in y or varies on scales much larger than R so that the Laplacian term 
remains indeed negligible. The zero order Lagrange parameters a', C, 7, involved in this 
expression, can be obtained by directly maximizing the entropy by the same method as in 



section (2.3.1). A relation between the jet curvature and topography is then obtained at 
first order. This approach is developed in appendix (|A|). 

However it is more simple to proceed differently : we start from the jet equation and 
show that its integrability condition provides the relation between the jet curvature and 
topography. To catch this effect, we take into account the radius of curvature of the jet, 
denoted r, like in section 2.3.4, but state that r is constant across the jet, assumed thin. 
From the Gibbs state equation (|56|), using the boundary layer approximation, we thus 



obtain the jet equation in term of the transverse coordinate ( : 

d 2 ip 1 dtp if) ( , ip 



de e ~ri + ^- Ky) = B - tanh r " c ^ + 1V ) ■ m 

We have introduced e = ±1 to account for the direction of curvature (keeping r > 0). We 
define e = 1 (respectively -1) if the curvature of the jet is such that 4>i ( respectively 1 ) 
is in the inner part of the jet. Note that, in the case of a vortex, as in our notations tp is 
proportional to the opposite of the pressure the case e = 1 ( resp e = — 1 ) corresponds to 
a cyclone ( resp an anticyclone ). 

The algebraic equation (^) depends on three Lagrange parameters, instead of two for 
the periodic case of previous section, but we have three additional constraints, the condition 
that ip has the same value at the two boundaries, the circulation constraint ( |52j ) and the 
momentum constraint (|54"|), This will be achieved in general by boundary jets. Let us first 
study the interior jets. 

To study the interior jet, we make the change of variables : 

t =r-^ = — — + m (61) 
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We assume the variations of y in the jet width are negligible ( R « scale of variation of 
h(y) ), so that y is treated as a constant. Then we obtain the jet equation : 

_ B. # + + a^w =B+ + 

with (j) — > (j)±i for t — > ±00, where again <p±i corresponds to the solutions of the algebraic 
equation (|59"D, rescaled as 

<t> + ^-B + ^-h{y)= tanh(C^) . (63) 



Let us consider, as in section ( p. 3.4 ), the analogy of the equation (|62|) with the motion 



equation of a particle in the potential : 



UW = _ £ + ( B + - 2l±2» ) . .. 



Like in section 2.3.4., integration of (|62j ) from —00 to +00 imposes the integrability condition 



-00 



The second term of the l.h.s. of equation (^) can be interpreted as a friction term: if e = 1, 
the 'particle' starting from rest at 4>\ can reach a state of rest at 0_i only if the difference 
of 'potential' corresponds to the energy loss ( p5| ) by friction (if e = —1 the same is true in 
the reversed direction). 

As in the periodic case we have made the thin jet assumption R « r, so that the 
friction term ( rhs of (|65|)) is a correction of order R/r : U(<fix) — U((f)-i) = O(Rfr). We 
first neglect it to get the zero order results, so we write U{4>) = U(4>—i). Therefore the 
two hatched areas in figure ^ must be equal, like in figure Due to the symmetry of 
the tanh function, this clearly implies that the central solution of the rescaled algebraic 
equation ( p3| ) must be 4>q = 0, so that a' /Co — B + 7y/Co — h(y) = (denoting the zero 
order Lagrange parameters by the index 0). This is possible at different latitudes y only if 
ly/Co — h{y) = 0, or is of order R (so that it can be neglected at zero order). Then the 
integrability condition becomes 

a' = C B. (66) 

Furthermore <p±i are symmetric with respect to 0, of the form <f>±\ = ±u, determined by 
equation (|33|) , like in section 2.3. This parameter u is again related to the energy by (|32]). 
Finally, the terms 72//C0 — h(y) = and the curvature term disappears in the jet equation 
(|62l) , which therefore reduces to (^), discussed in section 2.3.2j . 



The first order solution outside the jet is obtained as a small correction 5(f>±i to the 
zero order solutions iu, with also a small correction a% and C\ to the parameters Oq/Co 
and Cq, 



a' a' n 



f> ±1 = ±u + 5cj> ±1 {y) , ^ = ^ + ai , C = Co + d (67) 
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From (|65|), we deduce that U((pi) — U{<j>-\) has the sign of e. As may be seen in figure 
|2], when a\ is positive, the line <t> + jj moves upward, so that U{4>\) < U(4>-i). Thus a\ 
has the sign opposite to the sign of U{4>\) — U((j)-i) ; and we conclude that — ea± is always 
positive. Introducing this expansion ( |67| ) in the algebraic equation (|63|), using the zero 
order results (^) and (|66|), we obtain : 



r XM , s -aigp + d[B± Coup. - u 2 )] + C h(y) 

Co6J> ±1 (y) = 1-CM-ui) (68) 




Coming back to the stream function ip, using (61), we deduce the corresponding velocity 
v by differentiation with respect to y. This velocity outside the jet is zonal, along the unit 
vector x, and verifies: 



(69) 



It is therefore a constant plus a term proportional to the local beta-effect Notice that 
the corresponding shear ^~ is stronger than the deep shear d^/dy 2 = R 2 d 2 h/dy 2 by the 
factor [1 — Co(l — u 2 )]^ 1 > 1. The integrability condition (|65| ) now provides the curvature 
of the jet. We can approximate the r.h.s. of this relation, of order R, by the zero order jet 
profile (j36|), denoting : 



( see figure y(b) ). The l.h.s. can be expanded, using (|67|) . We first expand the ex- 
pression ( |64[ ) of the potential, U(<j>) = Uo(<p) + C\/Cq [^>tanh(Co0) — log cosh(Co«^)/Co] + 
[h(y) — jy/Co — ai]<f). We can approximate <j) ~ ±u in the correction terms, and ex- 
pand Uo(4>) = Uo(±u) + dUo / d(ft(±u) . The zero order equilibrium condition imposes that 
dUo/d(f>(±u) = 0, so that (p5|) becomes 

eu (h(y) - ^ - aA = e(u)^ (71) 

This equation (fn]), expresses the dependence in latitude y of the curvature radius r of the 
curve on which the jet is centered, thus defining the shape of the sub-domain interface as 
a function of the topography. 

Without topography and for 7 = 0, we get a constant jet curvature. The same result 
was obtained in sub-section 2.3.3 by a different argument of free energy minimization. The 
parameter u, related to the energy by (|3^) and to Co by (|33"|), quantifies the strength of 
the jet. By contrast, the vortex area is determined by the constraint on PV patch area 
(parameter B), but it is also related to the jet curvature, proportional by fl7l]) to the small 
shift ot,\ in chemical potential and temperature. Likewise the equilibrium temperature at a 
liquid-gas interface slightly depends on the bubble curvature, due to capillary effects. 



As explained in the end of section (3.2) the parameter 7 is linked to the zonal propaga- 



tion speed of the structure. The term 72/ in (|7l|), combined with a usual beta effect (linear 
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topography term h(y)), leads to an oscillation with latitude y of the jet curvature 1/r, i.e. 
a meandering jet. Another possibility is an exact compensation of the beta-effect by the jy 
term, leading to a propagating circular vortex, and the selection between these two alter- 
natives is discussed in next sub-section. An oval shaped, zonally elongated vortex, such as 
on Jupiter, is obtained when this compensation occurs, but with an additional quadratic 
topography h(y). Indeed, to get a zonally elongated vortex, supposed latitudinally centered 
in zero, the radius of curvature of the jet must decrease for y > and increase for y < 0. As 
a consequence, we deduce from ( |7l| ) that the topography must be extremal at the latitude 
on which the vortex is centered (it actually admits a maximum in the cyclonic case and a 
minimum in the anticyclonic case ). Moreover, we deduce from (^) that the surrounding 
flow must have a zero velocity at the latitude on which the vortex is centered and that 
the shear is cyclonic when the vortex is a cyclone and anticyclonic when the vortex is an 
anticyclone. More generally, the curvature can be related to the zonal velocity outside the 
jet, eliminating the topography between (|72]) and ([H]), 



Let us recall our approximations. Writing down the jet equation (|60|), when making 
the boundary layer approximation, we have assumed R « r. We also assumed that in 
the jet width the topography can be considered as a constant. If 1/y/a is a typical length 
scale for the topography variation this gives aR 2 « 1. Moreover we have assumed that 
the effective topography effect h{y) — jy/Co remains small along the jet. If Ly denotes the 
jet extension ( for example a vortex latitudinal size ) this approximation is valid as long as 

2 
V 



aLt « 1 



3.3.1 Beta effect or linear topography. 

Let h(y) = —(3y in the following of this section. (3 may mimic the beta-effect or a uniform 
velocity in the sublayer ( but we will refer it as the beta-effect ). 

A first class of equilibrium states corresponds to a single solution if)(y) of the algebraic 
equation (^). This determines a smooth zonal flow, with possibly intense jets at the 
boundaries y = ±g. The solution depends on the unknown parameters C, a' and 7, which 
are indirectly determined by the energy E, momentum M, and the condition (ift) = 0. The 
limit of small energy corresponds to C — > 00, for which we can neglect the term tp/R 2 on 



the left hand side of (59), which then reduces to if) = i? 2 /C[argtanh(/3y — B) + 72/ + a']. 
This corresponds indeed to arbitrarily small values of ip (small energy) as C — > 00. 

When the particular energy value E = R 2 (3 2 /2A is reached, a uniform PV is possible, 
with ip/R 2 = —j3y. Then PV mixing is complete, which clearly maximizes the mixing 
entropy. In this case, 7 = — Cf3, so that jy cancels the term Cip/R 2 in (p5S|). Physically, 



the uniform westward zonal velocity, 

v m = -R 2 f3 (73) 

tilts the free surface with uniform slope by the geostrophic balance, and the corresponding 
topographic beta-effect exactly balances the imposed beta-effect. 
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For a still higher energy, a first possibility is that again 



7 = -CP (74) 
so that the beta effect exactly balanced by the jy term in the jet equation of previous 



subsection. This cancellation is directly obtained in the general Gibbs state equations (55) 
and (|5"6|). Indeed the modified stream function ip' = ip + R 2 f3y satisfies the same equations 
as in the doubly-periodic case. Therefore in the limit of small R, the Gibbs states are made 
of subdomains with uniform ip' (uniform PV), separated by straight zonal jets or circular 
vortices. 

However jets or a vortex persist in this sea of uniform PV due to the constraint of energy 
conservation. The vortex moves westward at the same velocity v m , according to (|5^) , so 
they are just entrained by the background flow, without relative propagation (this can be 
physically understood by the cancellation of the beta-effect). 

The selection of the sub domain areas and PV values is given like in the periodic case 
of section 2, just replacing ip by ip' = ip — (3R 2 y. Therefore we get again probabilities 
p±i = (1 ± u) in the two subdomains with respective areas A±\ given by (|3l|), and stream 
function, 

ip ±1 = R 2 (B±u)- R 2 (3y (75) 

From this relation, we can calculate the energy E = i^iA-i + ipfA^/R 2 , so the energy 
condition ( |32| ) then becomes 

E = ^{u 2 -B 2 ) + ^- (76) 

Therefore these solutions with canceled beta effect can be obtained only beyond a minimum 
energy i? 2 /3 2 /24, corresponding to the potential energy of the surface tilting associated with 
the drift velocity v m . Then the excess energy will control the organization in two uniform 
PV areas. 

The shape of these subdomains can be obtained again by minimization of the jet free 
energy. However, unlike in the periodic case, jets occur at the boundaries y = ±^ as well 
as at subdomain interfaces. Indeed, such boundary jets are in general necessary to satisfy 
(v) = 0, or equivalently that the stream function ipi must be equal at the two boundaries 
y = ±i. In particular, the solutions (f75|) necessarily involve a stream function difference 
(or mass flux) —R 2 [3 associated with the drift velocity v m . This stream function difference 
must be compensated by boundary eastward jets with opposite total mass flux. We show in 
Appendix || that for two PV levels with similar initial areas a single eastward jet, separating 
two regions of uniform PV and weak westward drift, is the selected state (instead of two 
opposite jets in the periodic case). In the case of a strong PV level with a small initial area, 
the system organizes in a circular vortex like in the periodic case. In the limit u —> \B\, 
as one of the areas A±i goes to zero, the jet approximation falls down. The corresponding 
analysis of axisymmetric vortices and of the linear approximation for the Gibbs states, as 
performed in section 2, is still valid here. 

Up to now we have ignored the constraint of the momentum M (^J). This constraint 
imposes the latitude yo of the equilibrium structure ( a circular patch or a zonal band with 
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uniform PV ). For instance in the case B > 0, for which A\ > A—i (as seen from (pl|)), we 
define yo = f A yd 2 r/A\. Then 

M = I yqd 2 r = [ y(B + u)d 2 r + / y(B - u)d 2 r = 2uy A 1 (77) 

JD J A! J A-i 

We thus deduce the latitudinal position of the equilibrium structure : 

M , , 

* " 2uA7 (78) 

In the case of a single eastward jet, the subdomain position has been already fixed by 
the area (yo = A\/2). Then the only possibility to satisfy a moment M different from uA 2 
is that the jet oscillates in latitude with some amplitude A (then M — uA 2 ~ uA 2 ). This 
is possible if 7 7^ — (5C according to (|7l|), which becomes 

- = b(y- y ) (79) 
r 

where b = —(11(7 + 0C))/(2Ce(u)R) < and yo = a\jb. This equation clearly leads to a 
jet oscillating around the mean latitude yo ( as the curvature r is positive for y < yo and 
negative for y > yo ; recall that the curvature is by definition positive when positive PV 
is in the inner part of the jet ). Note that this oscillation propagates eastward at speed 
R 2 jj ( given by (|58|) ). Since b < 0, ^ > — /3, this speed is eastward with respect to the 
background drift v m (f7S|). 

3.3.2 Quadratic sublayer topography 



As explained in section ( |3.3| ), in the limit of small Rossby deformation radius R, the Gibbs 
state equation has solutions consisting of a vortex bounded by a strong jet on the scale of 
R. This corresponds to the case of an initial patch with strong PV and small area ( the 
asymmetry parameter B is sufficiently large ) with an energy sufficiently strong to get a 
structure of closed jet ( see figure ^). In the presence of a moderate topography h(y), this 
internal jet is no more circular but its radius of curvature r « R depends on y according 



to (|7T|). We have seen in previous subsection that a linear topography h(y) = —j3y leads to 
jets oscillating ( or to circular jets when 7 = —C/3 ). We shall discuss here how a quadratic 
term in h(y) modifies the shape of closed jets. We therefore assume a topography h(y) of 
the form : 

h(y) = ay 2 + by (80) 

This corresponds to a uniform deep zonal shear, with velocity Vd = R 2 d(h — [3y)/dy = 
2aR 2 y + b-f3. We focus our attention on vortex solutions, seeking close curves solutions of 
equation (|7l|). The vortices will be typically oval shaped as the ones seen on Jupiter. We 
then study how this shape ( for instance the ratio of the great axis of the oval to the small 
one ) depends on the topography ( sublayer flow ) and on the jet parameters. Application 
of these results to Great Red Spot observations will be discussed in next section. 

To make equation (|7l| ) more explicit, let s be a curvilinear parameterization of our 
curve, T(s) the tangent unit vector to the curve and 0(s) the angular function of the curve 
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defined by T(s) = (cos 0(s), sin 9 (s)) for any s. Then the radius of curvature r of the curve 
is linked to 0(s) by 1/r = d9/ds and ( |7l| ) yields the differential equations : 

= -dY 2 + 1 

81 

= sin 9(S) 

and ^ = C os0(S) (82) 

with c'X = x, c'Y = y — yo, c'S = s ; where d = (e 2 (u)R 2 a)/(eafu 2 ), yo = e(u)R(Cob — 
j) / (eatiCou) . The space coordinates X, Y and S are here non dimensional and have been 
obtained by dividing the real coordinates by the scale d = e(u)R/[eot\v). Note that as 
explained in section (|3.3|) ea\ > 0, so that d > 0. We further assume that a > 0, so that 
d > 0. 

We first note that the two variables and Y are independent of X. We will therefore 
consider the system formed by the two first differential equations (pl| ). It is easily verified 
that this system is Hamiltonian, with 9 and Y the two conjugate variables and 

# = cos - d— + y (83) 

the Hamiltonian. Thus is constant on the solution curves. We look for vortex solutions 
of our problem (81 and 82). Thus we require 9 to be a monotonic function of S. Moreover 
the curves must close, that is X and Y must be periodic. For symmetry reasons, it is easily 



verified that the solutions of (||,|8l]) with initial conditions 0(0) = f , Y(0) = ( H = ) 
and some X(0) are periodic. We prove in appendix (|C|) that these initial conditions are the 
only ones leading to closed curves. We also prove that the solutions of ( |8l| and ^) when 
d > d 

max = q does not define 9 as a monotonic function of S. They contain double points 
and thus are not possible solution for our problem. Once given these initial conditions, we 
can easily prove that the structure has both a zonal symmetry axis and a latitudinal one 
passing through yo- 

To study the shape of the jets, we numerically solve equations (^J8^) with initial 
conditions : 0(0) = f, y(0) = 0. We obtain closed curves with oval shapes, as shown in 



figure 1C. In figure O we have represented the width, the length and the aspect ratio of 



these vortices versus the parameter d. When d tends to g, the vortex width tends to a 
maximum value : u 
thus very elongated. 



9 ' 

maximum value : w max = | whereas the length diverges. In this limit, the vortices are 



4 Application to the Jupiter's Great Red Spot and Oval BC 

In previous sections, we have found maximum entropy states with the following properties 

• The fluid domain is partitioned in two subdomains with weak velocity, separated by 
jets whose width scales as the Rossby deformation radius. A strong initial PV level 
occupying a small area mixes in a subdomain with the form of a vortex bounded by 
an annular jet. 
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• In the presence of a parabolic topography h(y) (due to a sublayer zonal flow with 
uniform shear), outside the jet, exists a zonal flow (|69|) with uniform shear. The 
velocity at the latitude of the vortex center vanishes in the reference frame of the 
vortex. 

• The curvature of the jet is linked to the topography by (pfTj) . For the parabolic 
topography, solutions are oval shaped vortices, symmetric in latitude and longitude. 

These properties of our solutions are the main qualitative properties of the Jovian 
vortices. Moreover, would this description be correct, it would predict that the topography 
has an extremum at the center of the vortex. 

Dowling and Ingersoll (1989) derive the bottom topography, using the GRS and Oval 
BC velocity fields obtained from cloud motion. They analyze the results in the frame of a 
1-1/2 shallow water model (SW), with an active shallow layer floating on a much deeper 
layer. This deep layer is in steady zonal motion which acts like a topography h%. The 
SW topography h% is defined by fvd = —l/Ri(\)d(gli2)/d\ where Vd is the deep layer flow 
and Ri(X) is the latitudinal radius of curvature of Jupiter. Dowling and Ingersoll (1989) 
have deduced this SW topography hi by assuming the conservation of the shallow water 
potential vorticity + f)/h\ ( h\ is the upper layer thickness ) along the streamlines of the 
observed steady vortex flow. The vorticity uj is deduced from the measured velocity field, 
and the planetary vorticity / is known, so that the variation of h along each streamline is 
deduced. The pressure field is then obtained from the Bernoulli relation and the hydrostatic 
balance, leading to the field hi. The result depends on the radius of deformation of Rossby ( 
R* = gho/ fo, where ho is the mean upper layer height and fo the mean Coriolis parameter ), 
a free parameter in this analysis. Three test values have been chosen, R\ = 1700, R 2 = 2200 
and R% = 2600 km for the GRS and R\ = 1100, R* 2 = 1600 and R* 3 = 2000 km for the Oval 
BC ( we denote by a star superscript the physical parameters, to distinguish them from the 
non-dimensional quantities used earlier ). The height hi under each vortex has been found 
to depend only on latitude, and has been fitted as a quartic the planetographic latitude A: 

gh 2 = A + A x \ + A 2 X 2 + ,4 3 A 3 + A4A 4 . (84) 

The values obtained for the coefficients A{ in the vortex reference frame, for each of the 
vortices and for each of the values R^,R 2 , R3 are reported in table 1 of Dowling and Ingersoll 
(1989). 

Our model is the QG limit of this shallow water system. Starting from (SW) equations, 
we can derive the QG equations (|l]) by assuming the geostrophic balance and weak free 
surface deformation in comparison with the mean layer thickness ho. The validity of this 
QG approximation has been discussed by Dowling and Ingersoll (1989) and was found 
reasonably good as a first approach, although not accurate. We furthermore use the beta- 
plane approximation, linearizing the planetary vorticity around a reference latitude Ao ( 
Ao is taken to be —23° for the GRS and —33.5° for the Oval BC). Therefore we write 
f = fo + f3y, with fo = 20smAo and j3 = 2Vt cos \o/r z {\o) ( O is the planetary angular 
speed of rotation : 2ir/Vt = 9 h 55 mn 29.7 s and r z {\o) is the zonal planetary radius, which 
slightly depends on the latitude Ao, due to the ellipsoidal planetary shape, see formula (4) 
of Dowling and Ingersoll (1989) ) . We then obtain the QG potential vorticity (g) with the 
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QG topography h*(y*) linked to the SW topography (||) by: 



gh 



h *(y*) = Tjk+Py* (85) 



We have computed the QG topography d8^) using results of Dowling and Ingersoll 



(1989) for the SW topography (84) for the three values of the Rossby deformation radius 



R\,R\* and R%, for the GRS and for the Oval BC. The result in figure [12| shows that, for 
both the GRS and the Oval BC, the QG topography has an extremum at a latitude which 
is nearly the center of the vortex. As far as we know, this fact has not been noticed in the 
literature. This result is in agreement with the predictions of our model. We note moreover 
that the two extrema of the topography are minima, thus our model predicts anticyclonic 
shears around the GRS and the Oval BC, as observed. Figure |85] shows a comparison of 
the QG topography derived from Dowling results with a quadratic approximation, in the 
case R* = R2 = 2200 km. This shows that the quadratic approximation h*(y*) = a*y* 2 is 
a good approximation on the latitudinal extension of the GRS. This also provides values 
of the parameter a* © : a* = 9.2 10" 13 , a* = 7.2 10~ 13 and a* = 6.4 10~ 13 km" 2 s" 1 for 
R* = R\,R^ and #3 respectively. 

Let us deduce the corresponding non-dimensional parameters. First, the PV levels 
were normalized by (||), so our time unit T* will depend on the real PV level difference 
: {a\ — a*_i)/2 = y*. The other parameters are the Rossby deformation radius R*, the 
segregation parameter u and the topography coefficient a* (|80|). 

We will consider R* as a free parameter and use the following data from GRS observation 

• The jet width /*. Let us define the width of the jet I* as the width on which the 
jet velocity is greater than one half of the maximum velocity. We use velocity mea- 
surement within the GRS of Mitchell et al (1981). They have used small clouds as 
tracers to measure velocities, and have observed that the velocity is nearly tangential 
to ellipses. Using a grid of concentric ellipses of constant eccentricity, the velocities 
have been plotted with respect to the semi-major axe A of the ellipse on which the 
measurement point lies. Results have been fitted by a quartic in A and may be seen 



in figure 13. Using these results, assuming that the velocity profile in the jet is sym- 
metric with respect to its maximum, we choose as jet width I* = 5.6 10 3 km. In our 
model, the normalized jet width I = l*/R* can be computed from the jet equation 
([34|) as shown in figure ||(a). This determines the parameter u from the parameter 
R*. The corresponding theoretical jet velocity can be compared to observations in 



figure 13, for R* = 2500 km (the shape is not very sensitive to this parameter). The 



computed results for u versus R* is shown in figure 15 



The maximum jet velocity j)^. We will use the value ^ = 110 ms 1 ( Mitchell 
et al 1981). Using (|35|) and giving real dimension gives : 

. R* d<f>. . . 

V ma x = J^^l\max{u). (86) 
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d<p/dT\ max (u) has been obtained by solving the non-dimensional jet equation (|3 
and is shown in figure ^(b). This now determine T* from R*. 

The velocity shear surrounding the vortex. The ambient zonal shear measured at the 
latitude of the GRS from Limaye et al. (1986) is cr* = 1.5e — 5 s" 1 . Using (^) in its 
dimensional form, for a quadratic topography gives : 



2R* \ cosh 2 (C n) 



(87) 



This permits to compute a* as a function of R* (since u has been determined, as well 
as Cq, related to u by (|33|) ). 



The computed results for a* and T* versus R* are shown in figure [l4]. This shows that 
our determination of the topography is in agreement with the QG topography deduced 
from the shallow water model of Dowling and Ingersoll (1989) within a factor of two. The 
corresponding PV level difference a\ — a* 1 is comparable to the planetary vorticity fo at 
the latitude of the center of the GRS when R ~ 2400 km. For this value of R*, figure |l^ 
shows that u is very close to 1. Furthermore, as the GRS area is very small compared to 
the global area of a latitudinal band centered around the GRS, the non dimensional area 



occupied by the positive PV is very close to 1. Using this area expression (31) we conclude 
that B is very close to 1. Using the definition of B (JToj) , we conclude that a\ ~ and 
a* 1 ~ — /o- As discussed below a forcing mechanism by convective plumes incoming from 
the sublayer is expected to yield this result. 

The shape of the jet depends on the parameter d = ^g~2 a* -a* an< ^ on ^ ne l en §th scale 

c * = ^ . We can determine these two parameters from the previously determined values 
of a*,T* and u, and from the observed half width of the Great Red Spot : y^ax = 4900km. 
This permits to calculate c'*, a\ and d versus R*. Figure [0] show d versus R*. The dot line 
represents the critical value d = | beside which a vortex solution exists. The ratio of the 
length to the width of the GRS is approximatively 2, which would correspond to d = 0.441 
( figure ) ; this is very close to the critical value |. From figure |l6|, our model predicts 



that the Rossby deformation radius is R* = 1800 km. Figure 1C shows the actual shape of 
the vortex for d = 0.441. 

However, in the jet shape analyze, to obtain ( |3l| , |32^ ) we have supposed a*Ly << (a* — 
a* 1 )/2 where Ly is the maximal latitudinal extension of the vortex ( the topography part 
of the PV remains negligible with respect to the PV ). For the value of R* calculated above, 
we find a*y^ ax 2 = 3.910~ 5 s _1 whereas a* — a*_ l = 1.310~ 4 s _1 . We are thus at the limit of 
validity of our assumption. 

Now we can reverse the procedure and propose a predictive model of the Great Red 
Spot. Assume a steady deep zonal flow with uniform shear, = 2a*R* 2 (y* — y|), vanishing 
at an origin yj depending on our reference frame (which we shall choose in order to cancel 
the vortex drift). This flow is presumably generated by deep thermal convection but we 
are concerned here only with the dynamics of the upper layer, assumed stably stratified 
(due to cooling by radiative effects). We model this stratified upper layer as a shallow layer 



27 



with radius of deformation R* ~ 2400 km. This layer is submitted to a total beta effect or 
'topography' h(y) = ay 2 - (/3 + 2ay r )y. 

Assume that PV spots with value —fo, occupying small area proportion are randomly 
generated in this layer. This would be the result of intense incoming thermal plumes, as 
recently discussed by Ingersoll et al (2000), : conservation of the absolute angular momen- 
tum during the radial expansion leads to strong decrease of the local absolute vorticity, 
which comes close to zero. This means that in the planetary reference frame, a local vor- 
ticity patch with value —fo is created. The opposite vorticity is globally created by the 
subducting flow, but it is close to due the much larger area. This gives our time unit 
(2/0)- 1 and B = 1 - 2A. 

The outcome of random PV mixing with the constraint of the conservation laws is then 
a zonal velocity ( |69|) in the observed upper layer and a vortex with area ( |3l] ) with velocity 
profile shown by the dot curve in fig 14. The vortex moves with the upper velocity at y = 0, 
it drifts with respect to the deep layer at velocity (|73| ) so that the beta effect is suppressed. 
The shape is an oval symmetric in x and y, with aspect ratio computed from the shape 



parameter d ( figure 11 



Note that a slightly larger area , or stronger energy could lead to very elongated vortices. 
Then argument of free energy minimization show that this would lead to a single eastward 
zonal jet. This may explain the jet observed in the Northern hemisphere of Jupiter at the 
same latitude as the GRS. 

For smaller Jovian vortices such as the Oval BC, the size of the vortices is comparable 
with the Rossby deformation radius. Thus such vortices may be described as done for 
axisymmetric vortices ( section 2.3.4| ). This explain why such vortices do not have a 



quiescent core as the GRS. 

Let us describe dark brown cyclonic spots ( 'Barges' ) at 14 N on Jupiter. Their first 
interest is to stress that cyclonic vortices embedded in cyclonic shear exist on Jupiter. Let 
us go further. The greater of these barges is studied from Voyager observations by Hatzes 
at al (1981). The meridional velocities measured at the latitude of the center of the barge 
( Hatzes et al (1981), figure 7 ) show a boundary jet organization around the perimeter of 
the barge ( v max = 25 ms _1 ), see figure ||(b). The surrounding shear is such that the shear 
velocity at the maximum latitude of the barge is the same as the maximum jet velocity. 
Thus our approximation aLy « 1 is not good. We can however explain the elongated 



shape, similar to figure 1C (b) obtained for d very close to d n 



We conclude that the Gibbs state equation ( pq ) derived from maximization of entropy of 
the QG model (g) is in the limit of small Rossby deformation radius a model that explains 
the main qualitative features of Jovian vortices. The statistical mechanics itself explains 
the organization of a turbulent flow in coherent structures. 



5 Conclusion 

Our first result is to provide a general explanation for the emergence and robustness of 
intense jets in atmospheric or oceanic turbulent flows. In the absence of topography or 
beta-effect turbulence mixes potential vorticity in subdomains, and such jets occur at the 
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interface of these subdomains, with a width of the order of the deformation radius. From 
a thermodynamic point of view, this is like coexistence of two phases. Indeed the vortex 
interaction becomes short ranged in the limit of small deformation radius, and statistical 
mechanics leads to a thermodynamic equilibrium between two 'phases', with different con- 
centrations of the 2 Potential Vorticity levels. Another approach leading to the same result 
is to consider the general partial differential equation ( |l8|) characterizing the equilibrium 
states. This equation reduces to the algebraic equation ( po|) in the limit of small deforma- 
tion radius. The two uniform subdomains correspond to two solutions -0-1 and ip\ of this 
equation. At the interface of these subdomains, the general pde reduces to the equation 
(|36|), whose solution determines the jet profile. In addition, a solvability condition of this 
equation confirms the relation of equilibrium between the two 'phases', which was obtained 
in the thermodynamic approach. 

All our results have been obtained for a 2 potential vorticity level case, but cases with 
more levels would lead to qualitatively similar results, although the quantitative analysis 
would be more involved due to the additional parameters. 

In the presence of beta effect or topography, for low energy, purely zonal flows with 
gentle variation in latitude are obtained. A critical energy is the energy of the state where 
the zonal flow just compensates the beta-effect. For this state the PV is strictly uniform in 
the whole domain. For greater energy, two well mixed domain separated by jets appears, 
as in the without topography case. However the PV is no more strictly uniform in the well 
mixed subdomain : a zonal flow exists. Also the jet curvature depends on latitude. With 
ordinary beta effect this yields an intense eastward jet, purely zonal or wavy depending on 
the constraint on the momentum M. With the quadratic beta-effect generated by a deep 
shear, this can produce an oval shaped vortex. This vortex then drift so as to compensate 
the beta-effect. In other words, in the vortex reference frame the equivalent topography 
h{y) admits an extremum, and this is in agreement with the data of Dowling and Ingersoll 
(1989). 

Our quasi analytical approach therefore explains most of the basic features of the Great 
Red Spot and other Jovian vortices. It can be developed into a more accurate predictive 
model along the following lines. First the approximation R « r of thin jet is convenient 
for a qualitative understanding but is only marginally satisfied. However this limitation 
can be overcomed by numerical determination of the equilibrium state equation ( |i~8| ) by 
methods like used by Turkington and Whitaker (1995) or using relaxation equations toward 
the maximum entropy state as described by Robert et Sommeria (1992). Furthermore, 
extension to the more general shallow water model is desirable, as the Rossby number ( 
~ 0.36 where it a maximal ) is not very small. This can be formally achieved ( in preparation 

)• 

Finally the results rely on an assumption of ergodicity, or complete potential vorticity 
mixing consistent with the constraints on the conservation laws. Various numerical and 
laboratory experiments in the case of Euler equations ( see e.g. Brands et al. 1989 ) 
indicate that mixing may not be global but be restricted to active regions. Organization 
into local vortices, rather that at the scale of the whole domain, is more likely with a 
small radius of deformation, as vortex interactions leading to coalescence are then screened. 
This is observed for instance in the numerical computations of Kukharkin Orszag and 
Yakhot (1995). By contrast, the zonal shear here promote vortex encounters ( as observed 
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in Voyager data ) and we expect a much better relaxation toward the global statistical 
equilibrium, which involves always a single vortex in a given shear zone ( as it minimizes 
the interfacial free energy ). 
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A Determination of the Gibbs state by direct entropy max- 
imization, in the presence of a topography ( beta-effect ) 



In section (|3.3|), we studied the limit of small Rossby deformation radius in the Gibbs state 
equation ( |56| ) by considering the jet equation ( p0| ) and its integr ability condition (|65|). We 
deduced that the Gibbs states are composed of subdomains in which ip verifies the algebraic 



equation ( p9[ ) separated by an interfacial jet whose curvature verifies (71). The aim of this 



annex is to prove that these results can be obtained by directly maximizing the entropy, 



adapting the method used in section ( 2.3.1 ). 
Let us make the following assumptions : 

1. In the limit of small Rossby deformation radius, the probability p of finding the PV 
level ai takes two values p±i(y), depending only on y. We are looking for vortex 
solutions. The vortex shape is described by the length l{y) on which the probability 



p takes the value p~i(y) ( See figure [17 



The two subdomains where p take the two values p±\{y) are separated by a jet. The 
probabilities p±i(y) are supposed to be close to the their values without topography 
p±i = ±u, such that the free energy per unit length of the jet is well approximated 
by the one calculated without topography (|42|), If Ly denotes the vortex size, l/y/a 
a typical length on which topography varies, we will show this approximation to be 
valid as soon as aL v << 1. 

The boundary conditions can be relaxed ; that is no boundary term appears in the 
variation of the free energy at the order considered here. ( See discussion concerning 
boundary jets in section (|3.3|). 



Given these hypothesis, the Gibbs states is described by the 3 functions p±i(y) and 
l(y). We will determine them by maximizing the entropy S ( |l5| ) under the 3 constraints 
: energy d53), mass conservation ( |5ll) and momentum ([54]) . A necessary condition for a 
solution to this variational problem is the existence of 3 Lagrange parameters C,a and 7 
such that the first variations of the free energy : 

F=-S-^E + aM +1 M (88) 
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vanish. Using (15),(|53|),(|5l|) and (|||) ; @ and (|i~T|) and the above hypothesis ( see figure 



17 ), a direct calculation shows that the free energy (pq) is up to a constant : 



[f(pi(y),y)(l-l(y)) + f(p-i(y),y)l(y)]dy+ / j f( Pl (y),y)dy + 



f(pi(y),y)dy + LF Jet (u) 



Vvnax 



1\ 2 ,„„ , / 1 



with /(p,y) = (plogp+(l-p)log(l-p))-2C'o(p-~J -2(C £-a ) (p- - 
2(C %)-72/)fp-J 



(89) 



and where L is the jet length, Fj e t(u) is the jet free energy per unit length (42), calculated 
without topography. 

Considering first variations of the free energy ( |89] ) under variations of p\(y) ( resp p_i(y) ) 
proves that df /dp(p\(y),y) = and that df /dp(p-i(y),y) = 0. A direct calculation shows 
that : 

2 (p±i ~l)= tanh ( 2C o (p±i -5)+ C oMy) + CoB - o^j (90) 

Using ( pA[ ) and (^) ; recalling that we neglect the Laplacian term, a straight calculation 
shows that (0) is equivalent to the algebraic equation (|59|). 



Let us consider now first variations of the free energy (|89| ) under small variations Sl(y) 



of l(y). Using that the length of the jet is given by L = 2 f^™* Jl + \ (dl/dyfdy, a 
straightforward calculation shows that 5L = — fy max 5l(y)/rdy where r is the radius of 
curvature of the jet. We thus deduce from first variations of the free energy ( p9| ) : 

1 ' ,i{n) =f(pi(y),y)-f(p-i(y),y) (91) 



Hypothesis (2) : aLy « 1 permitted us to consider Fj et {u) as independent of y. In 
accordance with this hypothesis we evaluate f(pi(y),y) — f(p-i(y),y) at order zero, with 
P±i = ICli 1 ") at this order. We obtain : f(pi(y),y)-f(p-i(y),y) = 2u(al + C h(y) - jy) . 
Moreover using the free energy per unit length expression (]42]), ( |37|) and (|38|), one can show 



that Fj et {u) = 2e(u)CoR where e{u) is defined by (|70|). These two last results show that (91) 
is equivalent to ( [71] ) the expression for the radius of curvature r found by the integrability 
condition for the jet. 



B Boundary jets in the channel case. Beta-effect or linear 
sublayer topography. 

Let us derive the boundary jet properties in the case of a beta-effect (or linear sublayer 
topography). For the sake of simplicity, we will treat the case of zero circulation T = 0. 
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Let and be the values of <j> on the y = \ and V = —\ boundary respectively. The 
boundary jet satisfies the jet equation (|36|), but with boundary conditions <p(r = 0) = (j>^ 
and 4>{t — > +00) = 0i or Thus, using (|38|), we deduce that : 



\ (^r) 2 ( y = 4) + ^ = U{M = U{ ^ l] (92) 

(where the last equality comes from the integrability condition for the interfacial jet). 

We can relate dtfi/dr to the derivative dip/d£ normal to the boundary ( denoting the 
normal coordinate £ = ±5 =F V ), using (|6l|), 

f( y = 4)=L*4L = ± ±\^ a 1 (93 ) 



dT V 2 / Rd( \ 2 J C 
Furthermore the condition T = imposes 

,Ty dX + L,Ty dx = Tc{ y =2) + Tc{ y = -2 



2 



so that = \) = = -%). Then f|) becomes : 



2 



so that 



= W) (95) 



Note that the two values (/>^ cannot be equal. They must indeed satisfy, from ( |6l| ) and the 
condition of zero mass flux = if)— ), 

" C = -7/C (96) 

To solve these two equations ( |95|) and (|96|) , let us have a look at the potential U in 
figure 19. We have to compare — (p^ and 4>\ — 4>-\. We thus distinguish two cases. 
Using (|96|) and <j>±\ = ±u, we calculate : (4>^ — ^)/(</>i — 4>-\) = —j/(2Cu). We recall 
that u is of order unity. 

• Low beta-effect case : —j/C < 2u. 

For an x-independent statistical equilibrium, one possibility is a zonal band with 4>-i 
inside the band and near both boundaries, with <^ and (p^ at the boundaries. 
The symmetric solution, with 4^-1 an d 4>b-i a * the boundaries and 0_i near the 
boundaries, has the same free energy ( due to the symmetry of the potential U ). 
However the solution which maximizes the total free energy of the jets corresponds 
to 4>-l m the lower part of the domain, and 4>\ in the upper part, with 4to at the 
upper boundary ( y = \ ) and 4bo at the lower boundary ( y = — 5 ), and a single 
eastward interior jet. 
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High beta-effect case : j3 > 2u 

In this case, we note that — <jy^ > <pi — <j>-\. Then the two equations ( p6| ) and (95) 
determine a unique solution for 4>b) ( see n g ure 19 )• We are then necessarily in 
the case with (f> = <f>\ near the boundaries y = \ and <j> = (f>-\ near the other one. It 
again involves a single eastward internal jet. This jet can oscillate in latitude due to 
the momentum constraint ( according to (|79|) ). 



Once 4% are nx ed, we use equation (|6l]) on the two boundaries y = ±5 to conclude 
that a = — \{(f)fr + 4>b)- Then using the integrability condition ( |66| ) we calculate : 

<t = i?( B+ *±4) (97) 
This closes the determination of our parameters. 



C Boundary condition for the equation (|ST| , |82|) of the oval- 
shaped vortices boundary 



We want to see whether equations ( |81| , p2[) with d > 0, defining the curve formed by the 
jet, admit periodic solutions both in x and y corresponding to vortices, or not. As stressed 
in section ( |3,3.2 ), equations (]8l| ) derive from the Hamiltonian (|8^) ; y and 9 are the two 



conjugated variables. Let us study the phase portrait of H. For 9 in [0,2-7r[, there are 
4 critical points : Pi = (0,1/d), P 2 = (0,-1/d), P 3 = (tt, 1/d), P 4 = (tt,-1/(£). By 
linearization around these fixed points, one easily prove that P\ and P3 are stable fixed 
points whereas P2 and P4 are hyperbolic fixed points. This permits to draw the phase 



portrait : figure 18. Using H (^), we obtain that the unstable manifolds are given by 
1 — 2/(3v / d) = H and — l + 2/(3y / d) = H respectively. The parameter d governs a transition 
of the phase space structure. This transition occurs when the two unstable manifolds merge 
; this permits to compute the transition value for d : d = |. We are looking for vortices 
solutions of ( ^l|j82] ). We recall that 9 is the angle with the x axe. We thus impose 9 to be a 
monotonic function of s on a trajectory. Thus areas c of the phase portrait (figure |l8| ) are 
forbidden ; they would correspond to 9 varying in a finite interval. We won't address their 
analyze there, but oscillating jet solutions can be found in these areas. Area 62 : it can 
be shown that, as on such trajectories, as 9 is not strictly increasing, double points exists, 
thus forbidding area b2. Conversely areas a and bl are admissible trajectories, giving y as 
periodic functions of #.We have to had the condition that x ( |8"2| ) must also be a periodic 
function of 9 for the curve to close. Let us denote Ax the x variation when 9 runs in [0, 2tt]. 
We thus impose the condition Ax = 0. Using (p2]) and (pTl) we calculate : 



A f L m cos9de / f d[y 2 (9)-y^9-7r)} 

Ax = / cos 9ds = / — ; = / cos 9- — — - — — ^— r —d9 = U 

Jo Jo -dy 2 {9) + \ J_* (-d y 2(9) + l)(-dy 2 (9-ir) + l) 

(98) 

( The last expression is obtained rewriting the integral as a sum on [— f, f 1 plus a sum on 
^1 ; and performing a variable change ). Let us study the sign of y 2 {9) — y 2 {9 — it). 
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Using (H) we deduce that -d^+y{9) = H-cos6 and -d y ^ + y(9-n) = H + cos9. 
From these two relations we conclude that y{9 — n) = y(9) implies cos(#) = and that 
y{9 — 7r) = —y{9) implies H = 0. Thus if H / ; y 2 (9) — y 2 (9 — tt) does not change sign 
on [— f , §]■ Moreover, on the areas a) and bl) (—dy 2 (9) + 1) does not change sign. Thus 
if H 0, the argument of the last integral of (|98| ) does not change sign and Ax can not 
be zero. We thus conclude that the only solution where x is a periodic function of 9 is 
the solution corresponding to H = 0. This solution is the one obtained from (|8l|j8^ ) with 
initial conditions y(0) = and 9(0) = §• As we have previously excluded the area b2 ( when 
d > g ) we conclude that no vortex solution exists when d > | . We conclude that equations 
(81 ]82|) with d > 0, defining the curve formed by the jet, admit one periodic solution both 
in x and y corresponding to vortices, only when d < d max = |. This solution corresponds 
to H = (|83|). The vortex then admits a latitudinal and a zonal axes of symmetry. 
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Figure Captions. 



Figure || : Annular jets observed in the atmosphere of Jupiter, a) Velocity field in the 
Great Red Spot of Jupiter (20° South), from Dowling and Ingersoll (1989). b) Velocity 
field in the cyclonic Barge of Jupiter (14° North) from Hatzes et al (1981) 



Figure || : (a) Graphical representation of the algebraic equation fl20D, with the rescaled 
variable <fi = —a/C + ip/R 2 . The three solutions are at the intersection of the curve (left- 
hand side) and straight line (right-hand side). Here the integrability condition a = CqB 
for the differential equation 



Figure || : The free energy density f(p) ( |27| ) versus the probability p. For Co > 1 and 
(CqB — qo) small enough f(p) has two local minima and one local maximum, allowing to 
obtain two values p±% in the maximization of entropy under constraints. 

Figure ||| : The parameter u versus the Lagrange parameter Co, as the solution of (|33|), 

Figure ||| : Typical stream function profile in a jet (u = 0.75 ) versus the transverse 
coordinate r = C/R ( Top ) and corresponding velocity profile ( Bottom ) . 

Figure || : Jet properties versus the segregation parameter u. a) Jet width, defined 
as the width of the region with velocity greater than half the maximum jet velocity, b) 
Maximum velocity (d(j)/dT) max and jet kinetic energy e(u) ( dotted line). 

Figure [7| : Phase diagram of the Gibbs states versus the energy E and the asymmetry pa- 
rameter B. The outer line is the maximum energy achievable for a fixed B : E = ^-(1 — B 2 ). 
The frontiers line between the straight jets and the circular jets corresponds to A\ = 1/tt 
or A-i = 1/tt. it as been calculated using (||) and (||) : E = R 2 B 2 (2ir - 2)/(vr - 2) 2 . The 
dot line represents the frontiers between axisymmetric vortices and the circular jets. We 
define it as the energy value for which the circular vortex area A\ or A-\ ([}!]) is equal to 
(2l) 2 , where I is the typical jets width ( figure ^ ). Such a line depends on the numerical 
value of R the ratio of the Rossby deformation radius to the domain scale. It has been here 
numerically calculated for R = 0.03. 



Figure || : (a) Graphical representation of the algebraic equation (p0[), with the rescaled 
variable (f> = —a/C + ip/R 2 , like in figure ||, but in the case of a Gibbs state with an 
axisymmetric vortex (AC > 0). Then the rhs hatched area is greater than the lhs one. 
(b) The corresponding potential U((j)), given by (p9|), is asymmetric to compensate for the 
friction term in equation (|5C|), 

Figure |9] : Various axisymmetric stream-function profiles for decreasing AC ( AC = 
[0.90.60.30.10.050.030.01] and B = 0.75. 



Figure 1C : a ) Typical sub-domain shape with a topography h(y) = ay . The parame 



ter d has been chosen such that the ratio of the length on the width be 2 ; as on Jupiter's 
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GRS. b ) Typical sub-domain shape with a topography h(y) = ay when the parameter d id 
very close to its maximum value d = „. The shape is then very elongated, with latitudinal 
boundaries quasi parallel, as for instance the Jovian cyclonic vortices ( 'Barges' ) described 
by ( Hatzes et al 1981 ). 



Figure 11 : a ) Sub-domain non dimensional length and width versus the parameter d 
( topography h(y) = ay 2 ). b ) Sub-domain aspect ratio versus the parameter d. 

Figure |l2| : QG topography ( units ) versus latitude computed from data of Dowling 
and Ingersoll (1989) : a) under the GRS ; b) under the Oval BC. 



Figure 13 : Velocity profile within the GRS from ( Mitchell et al 1981 ). They have 
observed that the velocity is nearly tangential to ellipses. Using a grid of concentric ellipses 
of constant eccentricity, the velocities have been plotted with respect to the semi-major 
axis A of the ellipse on which the measurement point lies. Results have been fitted by a 
quartic in A. I* is the jet width, defined as the width on which the jet velocity is greater 
than half the maximum velocity. 

Figure 14 : A) Coefficient a for a quadratic topography h(y) = ay 2 versus R* com- 



puted from our QG model (|87|). The three cross show the coefficient a computed from 
QG topography deduced from Dowling SW observed results. B) Difference of PV lev- 
els a\ — a_i versus R* computed from our QG model (|86l) . The dot line represents the 
planetary vorticity Jq at the latitude of the center of the GRS. This show that one of the 
PV levels may be interpreted as vorticity generated by convection-plumes from the sublayer. 

Figure 15 : The segregation parameter u versus the Rossby deformation Radius R* for 
the GRS. u has been computed using the actual jet maximum velocity and width ( see 
section (f|) ). 



Figure 16 : The non dimensional parameter giving the shape of the curve : d ( see (81) 
) with respect to R* in our model of the GRS. The dot line represents the critical value 
d = | below which a vortex solution exists. The ratio of the length to the width of the 



GRS is approximately 2. From figure 11 we conclude that this correspond to d very close 
to the critical value |. From this figure, our model predicts that the Rossby deformation 
radius is R* = 1800 km ( see section (|j) for comments ). 

Figure [17] : Definition of l(y). 



Figure 18 : Phase portraits of the Hamiltonian H fl83| ) for yo = 0, governing the jet 
shape via differential equations (^|) ( two periods in 6 ). For vortices, we are looking 
for periodic solutions in y. Thus only trajectories of areas a) and b) are under interest. 
Conversely trajectories of area c) could correspond to oscillating jets. The parameters d 
governs a transition between two type of phase portraits. A) For d < | ( here d = 0.075 



trajectories of area a) can define y as a function of 9 corresponding to convex vortices. 

4 
9 



B ) For d > H ( here d = 0.075 ), for trajectories of area a), the curve y(9) admits double 
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points. Thus they can not define vortex boundaries. 



Figure 19 : Resolution of the equations (94 and |96] 
represent the two cases discussed in appendix (O). 



The long-doted and the doted lines 



39 



Figure 1: Annular jets observed in the atmosphere of Jupiter, a) Velocity field in the Great 
Red Spot of Jupiter (20° South), from Dowling and Ingersoll (1989). b) Velocity field in 
the cyclonic Barge of Jupiter (14° North) from Hatzes et al (1981) 
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Figure 2: (a) Graphical representation of the algebraic equation (|20|), with the rescaled 
variable <fi = —a/C + tp/R 2 . The three solutions are at the intersection of the curve (left- 
hand side) and straight line (right-hand side). Here the integrability condition a = CqB 
for the differential equation (|36|) is verified, so the two hatched areas are equal, 
b) The corresponding potential U((j)), given by (|37|), integral from to of the difference 
between the two curves (hatched area in (a)). 
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Figure 3: The free energy density f(p) (|27|) versus the probability p. For Co > 1 and 
(CqB — ao) small enough f(p) has two local minima and one local maximum, allowing to 
obtain two values p±\ in the maximization of entropy under constraints. 




42 



j;ure 5: Typical stream function profile in a jet (u = 0.75 ) versus the transverse coordinate 
= (,/R ( Top ) and corresponding velocity profile ( Bottom ) 
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Figure 6: Jet properties versus the segregation parameter u. a) Jet width, defined as the 
width of the region with velocity greater than half the maximum jet velocity, b) Maximum 
velocity (d(f>/dT) max and jet kinetic energy e(u) ( dotted line). 
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Figure 7: Phase diagram of the Gibbs states versus the energy E and the asymmetry 
parameter B. The outer line is the maximum energy achievable for a fixed B : E = ^-(1 — 
B 2 ). The frontiers line between the straight jets and the circular jets corresponds to A\ = 
1/vr or A_x = 1/vr. it as been calculated using @ and (||) : E = R 2 B 2 (2ir - 2) /(it - 2) 2 . 
The dot line represents the frontiers between axisymmetric vortices and the circular jets. 
We define it as the energy value for which the circular vortex area A\ or A_i ([}!]) is equal 
to (2l) 2 , where I is the typical jets width ( figure ^ ). Such a line depends on the numerical 
value of R the ratio of the Rossby deformation radius to the domain scale. It has been here 
numerically calculated for R = 0.03. 
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Figure 8: (a) Graphical representation of the algebraic equation (|20D, with the rescaled 
variable (j) — —a/C + t/j/R 2 , like in figure ||, but in the case of a Gibbs state with an 
axisymmetric vortex (AC > 0). Then the rhs hatched area is greater than the lhs one. 
(b) The corresponding potential U ((/)), given by (f49|), is asymmetric to compensate for the 
friction term in equation (pQ). 




4> + a/co-B 
tanh(Co<|)) 
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Figure 10: a ) Typical sub-domain shape with a topography h(y) = ay 2 . The parameter d 
has been chosen such that the ratio of the length on the width be 2 ; as on Jupiter's GRS. 
b ) Typical sub-domain shape with a topography h(y) = ay 2 when the parameter d id 
very close to its maximum value d = g. The shape is then very elongated, with latitudinal 
boundaries quasi parallel, as for instance the Jovian cyclonic vortices ( 'Barges' ) described 
by ( Hatzes et al 1981 ). 
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Figure 11: a ) Sub-domain non dimensional length and width versus the parameter d ( 
topography h(y) = ay 2 ). b ) Sub-domain aspect ratio versus the parameter d. 
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Figure 12: QG topography ( units s 1 ) versus latitude computed from data of Dowling 
and Ingersoll (1989) : a) under the GRS ; b) under the Oval BC. 



1.0 
0.7 
0.4 
0.1 
-0.2 
-0.5 



h(X) ( loV 1 ) 









GRS 



-30 -26 

R = 1700km 

R = 1700km 
h(X) ( loV 1 ) 



X ( degree ) 
-22 -18 -14 

R = 2200km R = 2600km 

a = 9.257E-13 km A -2s A -l 




X ( degree ) 



-38 -36 
R= 1100km 

R = 2000km 



-34 -32 
R = 1600km 



-30 



50 



Figure 13: Velocity profile within the GRS from ( Mitchell et al 1981 ). They have observed 
that the velocity is nearly tangential to ellipses. Using a grid of concentric ellipses of 
constant eccentricity, the velocities have been plotted with respect to the semi-major axis 
A of the ellipse on which the measurement point lies. Results have been fitted by a quartic 
in A. I* is the jet width, defined as the width on which the jet velocity is greater than half 
the maximum velocity. 
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Figure 14: A) Coefficient a for a quadratic topography h(y) = ay 2 versus R* computed 
from our QG model (p7|). The three cross show the coefficient a computed from QG 
topography deduced from Dowling SW observed results. B) Difference of PV levels a\ — a_i 
versus R* computed from our QG model (|86|). The dot line represents the planetary 
vorticity /o at the latitude of the center of the GRS. This show that one of the PV levels 
may be interpreted as vorticity generated by convection-plumes from the sublayer. 
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Figure 15: The segregation parameter u versus the Rossby deformation Radius R* for the 
GRS. u has been computed using the actual jet maximum velocity and width ( see section 
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Figure 16: The non dimensional parameter giving the shape of the curve : d ( see (^Tj) ) 
with respect to R* in our model of the GRS. The dot line represents the critical value d = | 
below which a vortex solution exists. The ratio of the length to the width of the GRS is 
approximately 2. From figure |ll| we conclude that this correspond to d very close to the 
critical value |. From this figure, our model predicts that the Rossby deformation radius 
is R* = 1800 km ( see section (|j) for comments ). 
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Figure 17: Definition of l(y). 




Figure 18: Phase portraits of the Hamiltonian H (^) for yo = 0, governing the jet shape 
via differential equations ( plj ) ( two periods in 9 ). For vortices, we are looking for periodic 
solutions in y. Thus only trajectories of areas a) and b) are under interest. Conversely 
trajectories of area c) could correspond to oscillating jets. The parameters d governs a 
transition between two type of phase portraits. A) For d < ^ ( here d = 0.075 ), trajectories 
of area a) can define y as a function of 9 corresponding to convex vortices. B ) For d > § 
( here d = 0.075 ), for trajectories of area a), the curve y(9) admits double points. Thus 
they can not define vortex boundaries. 
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Figure 19: Resolution of the equations (^4| and |96|). The long-doted and the doted lines 
represent the two cases discussed in appendix (Q). 
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